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We present molecular line observations, made with angular resolutions of 
■'20", toward the filamentary infrared dark cloud G34.43+0.24 using the APEX 



OO : [C0(3 -^ 2), i3C0(3 -^ 2), 01^0(3 -> 2) and CS(7 -^ 6) transitions], Nobeyama 

O : 45 m [CS(2 -^ 1), SiO(2 ^ 1), C34S(2 ^ 1), HC0+(1 -^ 0), Hi3C0+(l ^ 0) and 

^ : CH30H(2 -> 1) transitions] , and SEST [CS(2 -> 1) and 0^^0(2 -^ 1) transitions] 

in 

o, 

O . similar to that of the dust continuum emission observed with 11" resolution 

(Rathborne et al. 2005) showing a filamentary structure and four cores. The cores 
have local thermodynamic equilibrium masses ranging from 3.3 x 10^ - 1.5 x 10^ 
rS \ Mq and virial masses from 1.1x1 0^-1. 5x1 0^ Mq, molecular hydrogen densities 

5^ \ between 1.8 x 10"^ and 3.9 x 10^ cm^^, and column densities > 2.0 x 10^^ cm~^; 

values characteristics of massive star forming cores. The ^^C0(3 -^ 2) profile 
observed toward the most massive core reveals a blue profile indicating that the 
core is undergoing large-scale inward motion with an average infall velocity of 
1.3 km s~^ and a mass infall rate of 1.8 x 10^^ Mq yr~^. We report the discovery 
of a molecular outflow toward the northernmost core thought to be in a very 
early stage of evolution. We also detect the presence of high velocity gas toward 
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each of the other three cores, giving support to the hypothesis that the excess 

4.5 /im emission ("green fuzzies") detected toward these cores is due to shocked 
gas. The molecular outflows are massive and energetic, with masses ranging from 
25 - 80 Mq, momentum 2.3 - 6.9 x 10^ Mq km s~^, and kinetic energies 1.1 - 

3.6 X 10^ Mq km^ s~^; indicating that they are driven by luminous, high-mass 
young stellar objects. 

Subject headings: ISM: clouds — ISM: individual objects (G34.43+0.24) — ISM: 
jets and outflows — ISM: molecules — stars: formation 



1. Introduction 

Infrared dark clouds (IRDCs) are cold (<25 K), massive (~10^ - 10"^ M0), and dense 
(>10^ cm^^) molecular clouds with high column densities (~10^^-10^^ cm^^) seen as a dark 
silhouette against the bright mid-infrared background emission (Perault et al. 1996; Egan 
et al. 1998; Carey et al. 1998, 2000; Hennebelle et al. 2001; Simon et al. 2006a, 2006b). 
A wealth of recent observations show that IRDCs are the sites of high-mass stars and star 
cluster formation (Rathborne et al. 2006, 2007; Pillai et al. 2006; Jackson et al. 2008; 
Chambers et al. 2009). Dust continuum observations reveal that IRDCs harbor compact 
cold cores with typical masses of ~120 Mq and sizes <0.5 pc (Rathborne et al. 2006). 
The masses and sizes of the IRDC cores are similar to those of the massive hot cores (e.g., 
Garay & Lizano 1999), but IRDC cores are much colder. Recently, Chambers et al. (2009) 
investigated a sample of 190 cores found toward 38 IRDCs by Rathborne et al. (2006), and 
classified them as "active" if they are associated with enhanced 4.5 ;um emission, the so- 
called green fuzzies, and an embedded 24 /im source; and "quiescent" if they contain neither 
of these indicators. They found 37 active cores and 69 quiescent cores. Further, active cores 
are usually associated with bright (>1 Jy) methanol and water masers (Pillai et al 2006; 
Wang et al. 2006; Chambers et al. 2009), and have high bolometric luminosities, indicating 
that they are forming high-mass stars (M > 8 Mq) stars. The quiescent cores are the most 
likely candidates for high-mass starless cores. 

G34.43+0.24 is an IRDC with a filamentary morphology, extending by ~9' from north to 
south in equatorial projection (9.7 pc at the source distance of 3.7 kpc; Faiindez et al. 2004; 
Simon et al. 2006b). This IRDC is in the list of Rathborne et al. (2006) and Chambers et al. 
(2009). It is located roughly 11' north of the ultra-compact (UC) H II complex G34. 26-1-0. 15 
(Molinari et al. 1996). Molecular line observations toward this cloud were first reported 
by Miralles et al. (1994), who from ammonia observations with 1.5' resolution found an 
elongated structure in the N-S direction with a total mass of 1000 Mq. Near the center of 
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the filament lies the IRAS point source 18507+0121 which has a luminosity of 3.4 x 10^ Lq. 

Bronfman et al. (1996) detected toward IRAS 18507+0121 strong CS(2 -^ 1) emission 
with broad line wings, indicating it is associated with a dense massive star forming region. 
Ramesh et al. (1997) observed the IRAS source in the HC0+(1 -^ 0), Hi3C0+(l -^ 0), 
CS(2 — > 1) and C^'^S(2 — t- 1) molecular lines, with angular resolutions of ~16", modeling the 
observed line profiles as due to emission from a collapsing hot core that is hidden behind a 
cold and dense envelope. From observations of 3 mm continuum and H^'^CO^ line emissions, 
with ~5" angular resolution. Shepherd et al. (2004) detected two compact molecular cores, 
separated by 40", toward the center of the filament. Millimeter continuum observations 
revealed the presence of four dust cores (Faiindez et al. 2004; Garay et al. 2004; Rathborne 
et al. 2005), labeled by the latter authors as MMl, MM2, MM3, and MM4. The MMl dust 
core, which corresponds to the northern compact molecular core of Shepherd et al. (2004), 
contains a deeply embedded luminous object surrounded by several hundred solar masses of 
warm gas and dust. Based on the weak 6 cm continuum emission and the lack of detection 
at NIR wavelenghts. Shepherd et al. (2004) suggested that the embedded object appears 
to be a massive B2 protostar in an early stage of evolution. The MM2 dust core, which 
corresponds to the southern compact molecular core of Shepherd et al. (2004), is associated 
with IRAS 18507+0121; with a NIR cluster of young stars with a central BO. 5 star (Shepherd 
et al. 2004); with an UC H II region (Miralles et al. 1994; Molinari et al. 1998); with a 
variable H2O maser (Miralles et al. 1994); and with CII3OH maser emission (Szymczak 
et al. 2000). The MM3 dust core is located close to the northern edge of the filament, 
~3.5' north of MM2. Garay et al. (2004) suggested that MM3 corresponds to a massive 
and dense cold core that will eventually collapse to form a high-mass star. The MM4 core 
is located ~30" south of MM2 (Rathborne et al. 2005). Table [1] lists the measured 1.2 mm 
dust continuum emission properties as measured by Rathborne et al. (2006). Rathborne 
et al. (2005) and Chambers et al (2009) reported an excess 4.5 /zm emission toward MMl, 
MM3, and MM4 cores, which could be produced by either ionized gas and/or shocked gas 
and acts therefore as a signature of current star formation (MM2 is classified by the last 
authors as a "red" core, namely with bright 8 /im emission). They also showed that the 
four MM cores are associated with 24 fim point sources (emission from warm dust), which is 
an indicator of accretion onto an embedded protostar. Water maser emission, a well-known 
tracer of low-mass and high-mass star formation, has been detected toward all four MM 
cores (Wang et al. 2006; Chambers et al. 2009). Chambers et al. (2009) detected Class 
I CH3OH masers, another well-know signpost of star formation toward MMl, MM2, and 
MM3 cores. Rathborne et al. (2006) identified nine 1.2 mm continuum cores with masses 
ranging from 80 to 1300 Mq toward G34.43+0.24. 

Shepherd et al. (2007) observed the central region of the filament, encompassing 
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the MMl and MM2 cores, at millimeter wavelengths, in C0(1 -^ 0), ^^C0(1 -^ 0), and 
C^^0(1 — ^ 0) at ~ 4" angular resolution and at near infrared wavelengths with the Spitzer 
Space Telescope. They discovered five massive outfiows, three of them associated with the 
MM2 core and the remaining two associated with the MMl core. From the Spitzer data 
they identified 31 young stellar objects in the filament, with a combined mass of ~127 Mq, 
and additional 22 sources that they believe could be cluster members based on the presence 
of strong 24 /xm emission. Cortes et al. (2008) observed MMl and MM2 in polarized ther- 
mal dust emission at 3 mm and C0(1 — ^ 0) line emission with BIMA array. Their results 
suggest that there is a magnetic field othogonal to the axis formed by the MMl and MM2 
cores. Recently, Rathborne et al. (2008) carried out high angular submillimeter observations 
toward the MMl core using the Submillimeter Array. They determined, using continuum 
emission images, that MMl remains unresolved with a size of ~0.03 pc and a mass of 29 Mq. 
Their molecular line spectrum shows the presence of several complex molecules, suggesting 
that MMl is a hot core. They also found an extended ^^C0(3 — )■ 2) structure which may 
be evidence of a rotating envelope surrounding the central high-mass star, perpendicular to 
one of the outfiows that Shepherd et al. (2007) detected toward MMl. 

In this work we report observations of 11 molecular transitions using the APEX 12-m, 
SEST 15-m, and Nobeyama 45-m telescopes toward the G34.43+0.24 cloud. Since different 
molecular species trace different density and kinematics regimes, the capability of observ- 
ing multiple spectral lines is particularly advantageous to determine the morphology and 
kinematics of star forming cores. Due to the high abundance of CO, the C0(3 — ?■ 2) line 
is an ideal tracer of the low intensity, high-velocity wing emission. The ^^C0(3 -^ 2) line, 
being more optically thin than CO (3 — )■ 2), traces, in addition to the ambient gas, the low- 
velocity gas of the molecular outfiow. The SiO(2 — )■ 1) line is usually associated with the 
presence of outfiows since its abundance is highly enhanced by shocks. Methanol transi- 
tions are usually used as tracers of dense gas, temperature probes and tracers of shocks, 
since this molecule is evaporated from grain mantles. The HC0+(1 — ^ 0) and CS(2 — > 1) 
lines are high-density tracers and thus are good probes of the dense envelopes surrounding 
protostars. The 0^^0(3 ^ 2), 01^0(2 ^ 1), Hi3C0+(l -^ 0), CS(7 -^ 6) and C34S(2 -^ 1) 
molecular transitions probe the inner regions of the cores. The last two species are weakly 
infiuenced by infall and outfiow motions, and are used as reference to estimate the central 
ambient cloud velocity. 

The main goals of these multi-line observations are to determine the morphology and 
kinematics of the gas within the molecular cores in the G34.43+0.24 filamentary IRDC, in 
particular to investigate the presence of outfiowing gas, in order to increase our knowledge 
of the evolutionary state of the cores and to achieve a better understanding of the overall 
process of massive star formation along the filament. This is the first survey of several 
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molecular lines mapping portions and the whole cloud. It gives us, for the first time, a 
comprehensive view of the molecular gas distribution at the scale of parsecs. 



2. Observations 

The observations were made using the 12-m APE}4j telescope, the 15-m SEST telescope, 
both located in Chile, and the 45-m telescope of the Nobeyama Radio Observatory (NRO) 
located in Japan. The molecules, transitions, frequencies and observing parameters are listed 
in Table H 



2.1. APEX Telescope 

With the APEX telescope we observed the C0(3 -^ 2), 1^00(3 -^ 2), 0^*^0(3 -^ 2) and 
CS(7 -^ 6) lines. The observations were carried out during three epochs: 2005 July, 2005 
October, and 2006 May- July. The telescope beam size at the frequency of the observed 
lines ranged from 17 to 19" (FWHM). We used the heterodyne receiver (APEX-2a), which 
provided a channel separation of 122 kHz and a total bandwidth of 1 GHz with 8192 spectral 
channels. 

In CO (3 -^ 2) and ^^C0(3 —7-2), we mapped the emission within a 100" x 680" region 
which cover a spatial extent similar to that exhibited by the 1.2 mm dust continuum emission 
map (see Garay et al. 2004). The spacing between adjacent points on a regular grid was 
20" and the position-switching mode was used. The on-source integration time per map 
position was typically ~40 s in both molecules. System temperatures were in the range 110- 
340 K in CO (3 — >• 2) and 200-300 K in ^'^C0(3 —?> 2), resulting in an rms noise per channel, 
in antenna temperature, of typically 0.27 K and 0.37 K, respectively. In C^^0(3 — )► 2) we 
mapped the emission within a 80" x 100" zone, covering the central region of the filament also 
known as the IRAS 18507-1-0121 region. System temperatures were typically ~190 K. The 
on-source integration time per map position was ~50 s resulting an rms noise per channel 
of typically 0.29 K in antenna temperature. In CS(7 — > 6) we mapped the emission within a 
40" X 120" region. System temperatures were typically ~200 K. On-source integration time 
was ~80 s resulting an rms noise per channel of typically 0.18 K in antenna temperature. 

The molecular line data were reduced according to the standard procedure using CLASS, 
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a GILDAS working group software. 



2.2. NRO 45 m Telescope 

With the NRO 45 m telescope we observed the CS(2 -^ 1), SiO(2 -^ 1), 0^^8(2 -^ 1), 
HC0+(1 -^ 0), Hi3cO+(l ^ 0) and CH30H(2 ^ 1) lines toward the central region of the 
filament. The observations were undertaken during 1996 February-April. The telescope 
beam size at the frequencies of the observed lines is 15-18". We used an acousto-optical 
spectrometer with 37 kHz resolution covering a total bandwidth of 40 MHz. 

The emission was mapped in all lines within a 80" x 140" region. The spacing between 
adjacent points on a regular grid was 20", using the position-switching mode. System tem- 
peratures were in the range between 160 and 390 K. On-source integration times per map 
position ranged from 300 to 600 s resulting in rms noise of typically 0.08 K per channel in 
antenna temperature in all lines. First-order polynomials were fitted in order to estimate 
and remove the baselines. Part of this data were used by Ramesh, Bronfman & Deguchi 
(1997) to model the kinematics at the peak position of IRAS 18507+0121. 



2.3. SEST Telescope 

With the SEST telescope we observed the CS(2 -^ 1) and C^^0(2 -^ 1) fines toward the 
northern part of the filament (MM3 core). The telescope beam size is 24" at the frequency of 
the C^®0(2 — > 1) line and 51" at the frequency of the CS(2 — > 1) line. The observations were 
carried out in 2003 March. We used the high-resolution acousto-optical spectrometers, which 
provided a channel separation of 43 kHz and a total bandwidth of 43 MHz. We observed, 
with 30" spacing between adjacent points on a regular grid, 25 positions within a 150"xl50" 
region. The on-source integration times per map position were typically ~130 s for CS and 
~70 s for C^^O. System temperatures were in the range 230-270 K for CS and 400-540 K for 
C^^O, resulting in an rms noise, in antenna temperature, of typically 0.07 K and 0.18 K per 
channel, respectively. 
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3. Results 

3.1. Filament and Dense Cores 

Figure [T] shows tlie ^'^C0(3 -^ 2) spectra observed witli APEX along the whole filamen- 
tary structure detected in the 1.2 mm dust continuum observations of Garay et al. (2004; 
see their Figure 4). Emission in the ^^C0(3 -> 2) line was detected across the whole mapped 
region. 

Figure [2] shows contour maps of the velocity integrated ambient gas emission (veloc- 
ity range from 52.9 to 61.1 km s^^, determined by inspection of the line profiles) in the 
i3CO(3 ^ 2), CS(7 -> 6), CS(2 ^ 1), SiO(2 ^ 1), Hi3C0+(l ^ 0) and 0^^0(3 -^ 2) lines 
observed toward the central ~1.5' x3' region of the filament. The maps are overlaid on a grey 
scale image of the 1.2 mm dust continuum emission observed with SIMBA/SEST by Garay 
et al. (2004). We find that the morphology of the molecular line emission resembles that of 
the dust continuum emission, showing two components: a northern component associated 
with the MMl core and a southern component associated with the MM2 dust core. There 
are, however, some differences in the morphology of the emission in the different molecular 
lines. These are likely due to chemistry and/or optical depth effects and/or differences in the 
regimes of density and temperature traced by the different lines. In particular, the emission 
in the CS(7 -^ 6) line is less extended than in the other transitions and runs nearly parallel 
to the north-south direction. 

Figure (31 presents contour maps of the velocity integrated ambient gas emission (ve- 
locity range from 55.8 to 60.6 km s~^) in three methanol lines [CH30H(2o — )■ Iq) A+, 
CH30H(2_i -^ l_i) E, CH30H(2o -^ h) E] observed toward the central region of G34.43+0.24. 
The maps are overlaid on a gray scale image of the 1.2 mm dust continuum emission ob- 
served with SIMBA/SEST. The morphology of the emission in the three methanol lines is all 
similar. They are, however, different from the morphology exhibited by the dust continuum 
emission and the rest of the molecular lines. The methanol peaks do not coincide with the 
dust peaks, on the contrary they seem to be slightly anticorrelated. This could be due to 
depletion of methanol within the central regions of the dust cores. 

Figure m shows contour maps of velocity-integrated ambient gas emission (velocity range 
from 56.1 to 61.5 km s"^) in the 1^00(3 -^ 2), CS(2 -^ 1) and 0^^0(2 -^ 1) lines observed 
toward the MM3 core, overlaid on a gray scale image of the 1.2 mm dust continuum emission. 
Whereas the morphology of the emission in the ^^C0(3 -^ 2) line is similar to that of the 
dust continuum emission, that in the 0^*^0(2 — > 1) line is somewhat different, showing a 
peak that is displaced from the peak of the dust emission. 



Table |3] gives the observed parameters of the molecular cores. The angular sizes, de- 
termined from the contour maps, correspond to the observed semi-major and semi-minor 
axes. The line center velocity, line width and velocity-integrated brightness temperature 
were determined from Gaussian fits to the composite spectra of the ambient cloud emission 
within each core. The cores are assumed to be at a distance of 3.7 kpc, which corresponds 
to the kinematical distance of IRAS 18507-1-0121 reported by Faiindez et al. (2004) and the 
kinematical distance of the IRDC given by Simon et al. (2006b). 



3.2. Spectra Observed Toward the Dense Cores 

There are notable differences between the line profiles observed in different molecular 
species, for different positions. In what follows we describe the spectra observed toward 
three selected positions (0", 20"), (20", -20") and (40", 200") associated with, respectively, 
the MMl, MM2 and MM3 dust peaks. 

Figure |5] shows the spectra observed toward the (0", 20") position, which corresponds to 
the MMl core. The CO (3 — )■ 2) profile shows, in addition to considerable self-absorption, the 
presence of both blueshifted and redshifted high velocity emission. The strong CO (3 -> 2) 
wing emission at redshifted velocities is also clearly seen in the SiO(2 — > 1), CS(2 — > 1), and 
HC0+(1 -^ 0) spectra. 

Figure |6] shows the spectra observed toward the (20",— 20") position (MM2 core). The 
C0(3 — ?■ 2) profile shows considerable self-absorption, and the presence of high velocity emis- 
sion at both blueshifted and redshifted velocities. The ^^C0(3 — > 2) line profile exhibits a 
strong blueshifted feature, with a peak at the velocity of 55.8 km s~^, and a fiat redshifted 
shoulder. The C^^0(3 — > 2) line shows strong emission toward blueshifted velocities, with 
a peak at the velocity of 56.1 km s~^, and a skewed profile toward redshifted velocities. 
A similar profile is seen in the H13qo+(1 -^ 0) line. The CS(2 -^ 1) and HC0+(1 -^ 0) 
lines show double-peaked profiles; the former exhibits a strong blue peak and a weaker red 
peak whereas the later shows an opposite asymmetry, namely a strong redshifted peak and a 
weaker blueshifted peak. The CS(7 — )■ 6) and C^'^S(2 -^ 1) lines exhibit symmetrical profiles, 
with peak line center velocites, determined from Gaussian fits, of 58.1 ±0.1 km s^^ and 57.4 
± 0.2 km s~^, respectively. 

Ramesh et al. (1997) modeled the observed CS(2 -^ 1), C^^S{2 -^ 1), HC0+(1 -^ 0) 
and II^'^C0"'"(1 -^ 0) profiles toward MM2 as being produced by a collapsing hot core hidden 
behind a cold intervening screen or envelope. We have modeled here the line profiles of 
^^C0(3 — 7- 2), C^^0(3 — 7- 2) and CS(7 — )► 6) using the simple model of spectral-line profiles 
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from contracting clouds of Myers et al. (1996). We concluded that the MM2 core is under- 
going large-scale inward motions with a infall velocity of ~1.3 km s~^ (see section 14.1 ■2p . 

Figure [7] shows the spectra observed toward the (40", 200") position, which corre- 
sponds to the MM3 core. The CO (3 — > 2) spectrum shows, in addition to considerable 
self-absorption, strong blueshifted wing emission up to a velocity of 42.2 km s~^. The 
^^C0(3 -^ 2) spectrum exhibits a strong Gaussian line, with a center velocity of 58.5 ± 0.1 
km s~^, and blueshifted wing emission up to a velocity of 53.6 km s~^. These observations 
strongly indicate the presence of an energy source within the MM3 core. The C^^0(2 — >■ 1) 
and CS(2 -^ 1) spectra exhibit nearly Gaussian lines profiles with peak line center velocities 
of 58.9 ± 0.1 and 58.8 ± 0.1 km s~^, respectively. 



3.3. Molecular Outflows 

3.3.1. Central Region 

Figure E] shows position-velocity diagrams of the C0(3 -^ 2), CS(2 -^ 1), SiO(2 -f 1), 
and ^^C0(3 — ?■ 2) emission, toward the central region of G34. 43-1-0. 24, along two cuts in 
the north-south direction passing through the peak position of the MMl (left panel) and 
MM2 (right panel) cores. Strikingly seen in this Figure is the presence of high velocity gas 
associated with the MMl and MM2 cores, most likely due to molecular outflows emanating 
from both cores. Particularly notable is the detection of emission in the silicon monoxide 
line, which is a powerful tracer of the presence of strong shocks. The abundance of SiO 
molecules is highly enhanced, with respect to ambient abundances, in molecular outflows as 
a result of destruction of dust grains by shocks, giving rise to the injection into the gas phase 
of Si atoms and/or Si-bearing species. Once silicon is injected into the gas phase, chemical 
models based on ion-molecule reactions predict a large abundance of SiO molecules (Turner 
& Dalgarno 1977; Hartquist et al. 1980). 

In the cut passing through the peak of MMl, the most prominent outflow emission is 
seen redshifted with respect to the ambient gas velocity, extending up to velocities of 80.7 
km s-i in CO (3 -^ 2) and 70.2 km s'^ in CS(2 ^ 1) and SiO (2 -^ 1). Weak or no outflow 
emission is detected in the ^^C0(3 -^ 2) line, indicating low optical depths. In the cut 
passing through the peak of MM2, both blueshifted and redshifted high velocity emission 
are clearly seen emanating from the MM2 core. In CO (3 — ?■ 2) the blueshifted and redshifted 
emissions extend up to velocities of 42.0 and 74.2 km s^^, respectively. Also seen in this cut 
is wing emission emanating from the MMl core as well as from the MM4 core, suggesting 
that MM4 also contains an outflow. 
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Maps of the velocity integrated line wing emission, in the CO (3 -^ 2), SiO(2 — t- 1) and 
CS(2 -> 1) lines toward the central 1.5'x3' region of G34.43+0.24, including the MMl, MM2 
and MM4 dust cores, are shown in Figure IHl In SiO(2 -^ 1) and CS(2 — > 1), the range 
of velocity integration of the blueshifted emission (blue contours or solid lines) is 49.0 < 
visr < 52.9 km s~^, whereas that of the redshifted emission (red contours or dashed lines) is 
61.1 < visr < 65.2 km s~^. In C0(3 — )■ 2) the range of velocity integration of the blueshifted 
and redshifted emissions are, respectively, 42.0 < visr < 53.0 km s~^ and 62.0 < visr < 74.0 
km s~^. In Fig. IH] we identify two, possibly three, outflows associated with the dust cores 
within the region considered. The northernmost outflow (labeled 0-MMl) is associated with 
the MMl core. It is roughly oriented in the NE-SW direction, with the peaks of the redshifted 
and blueshifted emission being located, respectively, SW and NE of the 3 mm continuum 
source detected by Shepherd et al. (2004; marked with a star). There is, however, a large 
spatial overlap between the blueshifted and redshifted emissions. 

The outflow associated with the MM2 core, (labeled 0-MM2), shows a complex morphol- 
ogy with the blueshifted and redshifted emission exhibiting a considerable spatial overlap. 
The peak position of the redshifted and blueshifted emissions are located ~7" south of the 
UC HII region (marked with a triangle) and are coincident within the errors. This suggest 
that the symmetry axis of the outflow is close to the line of sight. 

Outflow-MM4 is associated with the MM4 dust core (marked with a square). The 
blueshifted and redshifted emissions exhibit a large overlap in their spatial distribution, 
suggesting that the flow is unresolved or that its axis of symmetry is close to the line of 
sight. 



3.3.2. Northern Region 

Figure [TO] shows position- velocity diagrams of the C0(3 — ?■ 2) and ^^C0(3 — > 2) emission 
toward the northern region of G34. 43-1-0. 24, along a line in the north-south direction passing 
through the peak position of the MM3 core. Clearly seen in the CO (3 -^ 2) cut is the presence 
of wing emission emanating from MM3, with velocities extending up to 42.2 km s~^ in the 
blueshifted side and 69.2 km s^^ in the redshifted side. In the ^'^C0(3 -^ 2) line the wing 
emission is barely seen. 

Figure [TT] shows maps of the integrated line wing emission, in the CO (3 — )■ 2) and 
^^C0(3 -7- 2) lines, toward the MM3 core. The range of velocity integration of the blueshifted 
emission (blue contours) is 41.2 < visr < 56.1 km s~^ for C0(3 -^ 2) and 52.6 < visr < 56.1 
km s~^ for ^'^C0(3 — )■ 2); whereas that of the redshifted emission (red contours) is 61.5 < 
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visr < 69.4 km s'^ for C0(3 -^ 2) and 61.5 < visr < 63.8 km s"^ for i3C0(3 -^ 2). In the 
CO (3 — > 2) map, the peak positions of the redshifted and blueshifted wing emissions are 
located, respectively, NE and SW of the peak position of the milhmeter dust core. The red 
and blue wing emissions exhibit however a considerable overlap in their spatial distribution. 



4. Discussion 
4.1. Molecular Cores 

4.1.1. Physical Parameters 

From the observations of the ^'^C0(3 — > 2) and C^^0(3 -^ 2) lines it is possible to esti- 
mate the column density and the mass of molecular gas within the central region (MMl and 
MM2 cores) as follows. Assuming that the C^^0(3 — > 2) line is emitted under local thermo- 
dynamic equilibrium (LTE) conditions and its emission is optically thin, then the column 
density is given by (e.g., Bourke et al. 1997) 

iV(C-O) ^ 8.084 X 10- (^e. + 0.88)exp(15.81/T..) ff 

^ ' 1 - exp(-15.81/Te,) J J ^ ' ^ ' 

where ris is the optical depth of the C^^0(3 -^ 2) line, Tex is the excitation temperature, 
and V is in km s~^. The LTE mass can be obtained from 

^lte\ ^ ,^, ( ^im \ AH2/Ci80]\ ( dV (Te, + 0.88)exp(15.81/Te,) [ [ , _ 
-M^) = '■''' [2?f2^) [-SJ^^J Wc) l-exp(-15.81/T..) J J ^^^^^^^ 

(2) 
where Hm is the mean molecular mass per H2 molecule, mu is the mass of a hydrogen atom, 
D is the source distance, [H2/C^^0] is the H2 to C^^O abundance ratio, and Q is in arcmin^. 
The opacities in the C^^0(3 — > 2) line are obtained from the ratio of the observed brightness 
temperatures in the C^^0(3 — ^ 2) and ^^C0(3 — > 2) lines as follows (see Bourke et al. 1997, 
for derivation) 

1 _ exp i-ns/r) _ [Jis{Te.) - Jls{T^,^)] T^b(-CO) 



1 - exp (-ns) [Jl3(Te.) - Jl3(Tbg)] T^b(Ci80) ' 

where Tbg is the background temperature, T^b is the main beam brightness temperature of 
the hne, J is given by (the subscripts "18" and "13" refer to the C^^O and ^^CO isotopes, 
respectively) 

JiT) = ^ I (4) 



-12- 



and, r = rig/ria is the optical depth ratio 



111 

^13 



13C0 



jkT^JhB.s + 1/3) expjE^^JkT,^) [i - exp{-hu,JkTe,)] 
{kTe^/hBi8 + 1/3) exp(E /fcTe^) [1 - exp{-hu^JkTe^)] 



(5) 



where B is the rotational constant of the molecule, E^ = hBJ{J + l) is the rotational energy 
at the lower state, and J is the rotational quantum number. Assuming an [^^CO/C^^O] 
abundance ratio of 7.6, we derive average optical depths in C^^0(3 — )■ 2) of 0.10 and 0.18 
toward the MMl and MM2 cores, respectively. Using an excitation temperature of 30 K and 
an [H2/C^^0] abundance ratio of 3.8 x 10^ (Wilson & Rood 1994) we derive total masses of 
330 and 1460 M© for the MMl and MM2 cores, respectively. The derived core parameters 
are listed in Table HI 

Masses can also be estimated assuming that the clouds are in virial equilibrium. Ne- 
glecting magnetic fields and external forces, the virial mass, Myij., for a spherical cloud of 
radius R is given by (MacLaren, Richardson & Wolfendale 1988) 

My,r = B(^) (y^X Me, (6) 

\pc/ \km s ^/ 

where Av is the average line width in km s~^, R in pc, and i? is a constant which depends on 
the density profile of the cloud. Assuming that the cores have uniform densities {B = 210) we 
derive virial masses of 1.1 x 10^, 1.5 x 10^ and 1.4 x 10^ M© for, respectively, the MMl, MM2 
and MM3 cores. These values correspond to the geometric mean of the masses determined 
from the transitions given in Table [5l The molecular hydrogen densities were computed from 
the mass assuming that the cloud has a spherical morphology and a uniform density. 

Table [6] summarizes the core masses determined from different methods and authors. 
Cols. (2) and (3) give the LTE and the virial masses, respectively, determined in this work; 
col. (4) gives the LTE mass from Shepherd et al. (2007); and, cols. (5) and (6) give the 
masses derived from dust continuum observations by Garay et al. (2004) and Rathborne et 
al. (2006), respectively. In our estimate, if the excitation temperature is 25 or 35 K, the LTE 
mass is decreased or increased by ~20%, respectively. Given the uncertainties associated 
with each method - such as in the abundance, excitation temperature, dust temperature, 
and dust opacity - the agreement is good between the masses determined in this work and 
the masses determined from dust emission. Differences in the mass obtained by different 
authors using dust emission are mostly due to different dust temperatures and core sizes 
used in their calculations. There are, however, discrepancies with the masses determined by 
Shepherd et al. (2004) from Hi3co+(l -^ 0) observations, of 4000 and 5000 Mq for MMl 
and MM2, respectively, and a total cloud mass of 50000 Mq. We note that Rathborne et al. 
(2005) estimated from dust emission observations a total cloud mass of 7500 M© assuming 
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a dust temperature of 30 K. A possible explanation of the larger mass estimate of Shepherd 
et al. (2004) is that the actual abundance of HCO"'" is probably larger than the adopted 
one, due to an enhancement produced by the molecular outflows. Moreover, Shepherd et al. 
(2007) determined, from observations of the C^^0(1 — > 0) line using the OVRO array, masses 
of 75 and 690 M0 for MMl and MM2, respectively. These values probably correspond to 
the inner mass of the cores. We note that masses derived from molecular lines might not be 
tracing all of the emission, and their emission may be optically thick, specially C^^0(1 — > 0). 
However, the assumption that the C^^0(3 -^ 2) is optically thin seems to be adequate and it 
is supported for the inspection of the line profiles and values obtained for the optical depths. 

Optical depths and column densities toward the peak position of the MMl and MM2 
cores are given in Table [71 The C^^S(2 -^ 1) and CS(2 -^ 1) values were obtained in an 
analogous way to that described above for the C^^0(3 — )► 2) and ^^C0(3 -^ 2) lines, adopting 
a [CS/C^^'S] abundance ratio of 22.5 and an [H2/CS] ratio of 1 x 10^ (e.g., van der Tak et 
al. 2000). 



4.1.2. MM2: a Collapsing Core 

The ^^C0(3 -^ 2) line emission at the peak position of the MM2 core shows a bright 
blueshifted peak and a fiat redshifted shoulder. Such line shape can be modeled as due to 
inward motions (Mardones et al. 1997). The blue peak is produced by high excitation gas 
on the far side of the envelope, while the redshifted peak results from low excitation gas in 
the near side of the envelope. The self-absorption due to low excitation, low velocity gas in 
the near side of the envelope gives rise to the central dip (Evans 1999). An estimate of the 
infall velocity, Vi„, can be obtained using the two-layer infall model of Myers et al. (1996). 

Myers et al. (1996) made a simple analytic model of radiative transfer in order to inspect 
if the observed spectra show any evidence of infall motions and to provide an estimate of 
the characteristic Vm- They assumed two uniform regions of equal temperature and velocity 
dispersion a, whose density and velocity are attenuation-weighted means over the front and 
rear halves of a centrally condensed, contracting cloud with n/rimax = {r/rmin)~^^'^ and 
V/Vmax = (r/r^i„)^^/2. We modeled the profile of the ^^00(3 -^ 2), 0^^0(3 -^ 2) and 
CS(7 — > 6) lines toward the MM2 core using their prescription and obtained the following 
model parameters: the peak optical depth, tq; the kinetic temperature, T^; the infall velocity, 
Vin] and the non-thermal velocity dispersion, cTj^j,. Figure [12] shows the observed and the 
best model line profiles, indicating motions consistent with a "fast" infall, Vm ~ cr. The 
values derived by the model are shown in the top right corner in the same figure. We note 
that some of the parameters for each fit are different, for example, the infall velocities are 



-14- 



1.5, 1.3, and 1.0 km s~^ for ^^CO, C^^O, and CS, respectively. This decreasing order is in 
accord with the level of asymmetry in the lines. The most asymmetric line, a single peak in 
the blue and a red shoulder, has the largest values of l^„ (Myers et al. 1996; Fuller et al. 
2005). Although the velocity dispersion in the CS line is a ~ 2Vi„, the other two lines have 
a ~ Vin indicating a fast infall. The derived systemic velocities show some scatter, most 
likely due to the influence of the molecular outflow, which is more intense toward redshifted 
velocities. If the contribution of the outflow is removed, the blue peak becomes stronger 
than the red shoulder, and this will have the effect of increasing the derived infall velocity. 

Ramesh et al. (1997) modeled the profile of the CS(2 ^ 1), C34S(2 ^ 1), HC0+(1 -^ 0) 
and H^^C0^(1 -^ 0) lines toward the MM2 core using both the large velocity gradient (LVG) 
method and the simple model of Myers et al. (1996). For the LVG method they adopted a 
constant temperature of 22 K (we determined a kinetic temperature of ~36 K), and found 
a collapse velocity that varies radially as 8(r/ro)~'^'^ km s"^ with r^ = 4.5 x 10^^ cm. At a 
radius of r = 5.5 x 10^^ cm (equivalent to an angular size of 10") we get Vm = 1.8 km s^^. 
Given the uncertainties associated with each method, the agreement is good with the values 
determined by us. The infall speed estimated by Ramesh et al. (1997) using the Myers et 
al. (1996) analytic model is about 1 order of magnitude smaller than the velocities they 
determine from the LVG analysis; however, since the fit of the Myers model is admittedly 
not too good, they argue in favor of the larger value, which agrees well with ours here. 

The mass infall rate can be estimated assuming M = M^Vin/ R, where M/j is the mass 
within the radius R and Vin is the infall speed at radius R. Using the LTE mass derived 
at the central position of MM2, of 240 M0, an average infall velocity of 1.3 km s^^, and 
a radius of 5.5 x 10^^ cm (corresponding to a 10"angular radius); we obtain a mass infall 
rate of 1.8 x 10^^ Mq yr^^. Such accretion rate is high enough to overcome the radiation 
pressure of the central star and to form a massive star. Shepherd et al. (2007) identified 
two likely massive stars (M > 8 Mq) with the MM2 core, with one of them, located near 
the projected center of the UC H II region, having a spectral type of BO. 5 (Shepherd et al. 
(2004). Models of massive envelopes accreting onto this kind of young stars require accretion 
rates of M > 5 x 10~^ Mq yr~^ (Osorio, Lizano & D'Alessio 1999), in agreement with our 
results. 



4.2. Molecular Outflows 

To compute physical parameters of the outflowing gas, we followed the standard LTE 
formalism described in Bourke et al. (1997) using the emission in the C0(3 -^ 2) line. 
The main source of error in determining the mass in the outflow arises from the difficulty 
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in separating the contribution to the outflows of the emission in the velocity range of the 
ambient cloud. We adopt as velocity boundary between the blue and red wing emission and 
the ambient emission the values of, respectively, 53.0 and 60.8 km s~^ for the MMl outflow, 
51.8 and 60.9 km s'^ for the MM2 outflow, and 56.1 and 61.5 km s'^ for the MM3 outflow. 
Another source of error is the possibility that the CO (3 -^ 2) line might not be optically 
thin, particularly at low flow velocities. Fortunately, with the available data, we can asses 
this effect and make corrections. 

The mass in the high velocity wings is estimated from the C0(3— )'2) observations as 
follows. Assuming that the energy levels of CO are populated according to LTE and that 
the outflow emission is optically thin, then the total CO column density, A^(CO), at each 
observed position is given by (e.g.. Garden et al. 1991) 

AT^rn^ 7 7 ini3 (^e. + 0.92) exp(16.60/Te,.) 1 /"^ r W (7^ 

iV CO = 7.7 xlO — ^ ififin/-r M — TTFr~\ TFFT] / TB{v)dv , 7 

[1 - exp(-16.60/re^)J [J(Te^) - J{Tbg)\ J 

where Tex is the excitation temperature of the population in the rotational levels, Tbg is 
the background temperature, Tg is the brightness temperature of the C0(3— >-2) emission at 
velocity v, v is measured in km s~^, and J is defined by equation (j3]). The mass is then 
computed from the derived column densities as, 

M = [H2/CO] /i„ A ^Ar(CO), (8) 

where /i^ is the mean molecular mass per H2 molecule, [H2/CO] is the molecular hydrogen 
to carbon monoxide abundance ratio, A is the size of the emitting area in an individual 
position, and the sum is over all the observed positions inside each lobe. 

From the observations of the C0(3 — )■ 2) and ^'^C0(3 — )■ 2) lines it is possible to estimate 
the opacities of the flowing gas and therefore to assess whether or not the assumption that 
the CO (3 -^ 2) wing emission is optically thin is correct. Using expressions ([3]) and ([5]) for 
i^CO and ^=^C0, assuming T^^ = 30 K, and [^^co/^^co] = 50, we find that the average 
C0(3 -^ 2) opacities of the blueshifted gas (in the velocity range from 49.5 to 53.0 km s^^) 
and redshifted gas (in the velocity range from 60.8 to 68.9 km s~^) of the MMl flow are 
3.4 and 2.5, respectively. The given ranges are those in which emission is detected in both 
lines. For the MM2 outflow we find that the average optical depth of the blueshifted gas 
(in the velocity range from 47.2 to 51.8 km s~^) and redshifted gas (in the velocity range 
from 60.9 to 68.7 km s~^) are 5.4 and 4.4, respectively. For the MM3 outflow we find that 
the average optical depth of the blueshifted gas (velocity range from 52.8 to 56.1 km s^^) 
and redshifted gas (velocity range from 61.5 to 65.2 km s^^) are 7.1 and 4.2, respectively. 
These calculations indicate that the CO (3 — )■ 2) emission from the flowing gas, in the given 
velocity ranges, is moderately optically thick for the three flows. On the other hand, as the 
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position-velocity diagrams suggested, the ^^CO flow emission is optically thin, with optical 
depths between 0.05 and 0.13. 

To correct the mass determined under the optically thin assumption by the opacity 
effect we divided the CO (3 — > 2) wing emission into two regimes, an optically thin regime, 
in which no ^^C0(3 — > 2) emission is detected, and an optically thick regime. We assumed 
an [H2/CO] abundance ratio of 10"^ (Frerking, Langer & Wilson 1982) and an excitation 
temperature of the outflowing gas of 30 K. We derived masses in the optically thin regime 
of the 0.5, 1.1, and 1.6 Mq for the blueshifted lobes of MMl (38.6 < vi^r < 49.5 km s"^), 
MM2 (35.9 < wisr < 47.2 km s'^), and MM3 (40.1 < wisr < 52.8 km s"^), respectively, and 
masses of 0.7, 0.8, and 0.7 Mq for the redshifted lobes of MMl (68.9 < -wisr < 74.5 km s"^), 
MM2 (68.7 < -yisr < 74.5 km s~i), and MM3 (65.2 < -yi^r < 74.3 km s"^), respectively. 
In the optically thick regime the mass determined under the optically thin assumption is 
multiplied by the factor t^/{1 — e^'^"'), where r^ is the average C0(3 — ?■ 2) optical depth 
of wing emission in the optically thick velocity range. We derive masses in the optically 
thick regime, corrected by opacity, of 2.6, 13.6, and 8.5 Mq for the blueshifted lobes of MMl 
(49.5 < fisr < 53.0 km s"^), MM2 (47.2 < tJi^^ < 51.8 km s'^), and MM3 (52.8 < t;isr < 56.1 
km s~^), respectively. Similarly, for the redshifted lobes we derived masses of 11.9, 30.8, and 
7.3 Mq toward MMl (60.8 < v^, < 68.9 km s'^), MM2 (60.9 < v^,, < 68.7 km s'^), and 
MM3 (61.5 < t^isr < 65.2 km s'^), respectively. 

To estimate the contribution to the flow mass from the low velocity material emitting in 
the same velocity range as the ambient cloud gas, we followed the prescription of Margulis 
& Lada (1985). Using expression (A16) of Bourke et al. (1997), 

/ Tinb dvdVt= I { J {vr - Vb) dn , (9) 

lobe Jvi, Jlohe \ ■^ J 

where T'', v^ and T'', v^ are the brightness temperature and velocity at the blue and red 
velocity boundaries, respectively. We estimated that the hidden flow masses corrected by 
opacity in the low velocity range within the blueshifted lobes toward MMl, MM2 and MM3 
are 3.0, 20.6 and 13.9 M©, and within the redshifted lobes are 6.2, 11.9 and 8.9 Mq, respec- 
tively (see Table ED. Finally, the total mass of the outflows are: 24.9 Mq for MMl, 78.8 M© 
for MM2, and 40.9 M© for MM3. 

Using the standard LTE formalism, it is also possible to estimate the momentum, P, 
and the kinetic energy, i?k, in the flow as follows (Margulis & Lada 1985; Bourke et al. 1997) 



P = Iff M{v,a,6)vdvdad6 , 
Ek = - M{v,a,6)v'^dvdad6 . (10) 
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These equations assume no correction for flow inclination, and so are strict lower limits. 
Another method to compute the physical parameters of the outflowing gas is to assume 
that all the mass is flowing at a velocity characteristic of the entire flow, Khar- The flow 
parameters are, in this approach, given by 

P = M Khar , E^=]^M KLr , (H) 

where M is the total mass of the outflow lobe. The blueshifted and redshifted emissions 
from the outflows show considerable overlap in their spatial distribution, indicating that 
the symmetry axis of the flows is likely to be near the line of sight. Hence, we adopt as 
a characteristic flow velocity the average of the maximum observed outflow velocities. The 
parameters determined under this assumption are likely to correspond to upper limits. The 
maximum velocities observed toward the blue and red lobes are, respectively, 12.0 and 23.1 
km s"^ for the MMl outflow, 15.0 and 15.9 km s'^ for the MM2 outflow, and 17.3 and 
11.6 km s~^ for the MM3 outflow. Table [H] gives the derived parameters of the blueshifted 
gas and redshifted gas of the flows calculated using both methods. The geometric mean 
values for the outflow parameters are Pi ~ 230, P2 ~ 690, and -P3 ~ 270 M© km s~^; and 
Ek, ~ 1290, Ek2 ~ 3630, and Ekg ~ 1120 Mq km^ s^^ (subscripts 1, 2, and 3 refer to 
the MM cores). If the excitation temperature is increased to 50 K, the flow parameters 
increase by 7%. The derived values of the outflow parameters are similar to those of massive 
and energetic molecular outflows driven by high-mass young stellar objects (Shepherd & 
Churchwell 1996; Beuther et al. 2002; Wu et al. 2004). We do not provide estimates of the 
dynamical timescale, mechanical luminosity, and mass outflow rate of the outflows toward 
G34.43+0.24 because they are likely to be aligned close to the line of sight, and therefore 
the lengths of the lobes are difficult to estimate. 

Rathborne et al. (2005) reported an excess 4.5 /xm emission toward the MM cores. The 
Spitzer IRAC 4.5 yum band contains the H I Bra line (at 4.052 /im), the C0(1 ^ 0) P(8) 
band head (4.6-4.8yum), and many H2 vibrational transitions (the most important being 
H2(0 —7- 0) S(9) line at 4.694 /im), but unlike the other three bands it does not contain any 
strong polycyclic aromatic hydrocarbon (PAH) features (e.g.. Smith et al. 2006). It has 
been postulated that the main contributers of the excess 4.5 /im emission could be emission 
from the C0(1 — > 0) P(8) band head (e.g., Marston et al. 2004) and/or emission from the 
H2(0 — > 0) S(9) line (e.g., Noriega-Crespo et al. 2004). Thus, the enhancement of the 4.5 /im 
emission could be a probe of jets and molecular outflows (e.g., Teixeira et al. 2008; Smith 
et al. 2006; Noriega-Crespo et al. 2004); therefore of current star formation. Rathborne et 
al. (2005) first found evidence for the presence of high velocity gas toward the MMl and 
MM3 cores from the broad line widths (AV" ~10 km s~^) observed in the CS(3 — )■ 2) and 
HCN(3 -^ 2) lines. A detailed study of the central region of G34.43+0.24 was carried out by 
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Shepherd et al. (2007), who reported two molecular outflows emanating from MMl (G34.4 
MM, in their paper), and three outflows originating from the MM2 core. They estimate a 
total mass, momentum and kinetic energy of 34.8 Mq, 350 M© km s~^, and 4.9 x 10^^ ergs 
(2460 M0 km2 s'^), for the outflows in MMl; and 111.2 Mq , 830 Mq km s-\ and 7.8 x 10^'^ 
ergs (3920 Mq km^ s~^), for the outflows in MM2. The masses are ~20% greater than our 
estimate, and the momenta and kinetic energies are within our estimated range. In spite of 
the different molecular transitions, angular resolutions and assumptions used. Shepherd et 
al. and our results are in agreement, within the observational errors. 



4.2.1. The Discovery of a New Molecular Outflow Toward the Core 034-43+0.24 MM3 

We report the discovery of one massive molecular outflow associated with the MM3 
core discovered by Garay et al. (2004). They suggested that this massive and dense core, 
without counterpart at mid-infrared (MSX) and far-infrared (IRAS) wavelengths implying 
that it has a low temperature, seemed to be in a quiescent state, and thus appeared to be 
the youngest of the massive cores in the filament. The presence of 24 /im emission, a green 
fuzzy, water and methanol masers (Rathborne et al. 2005; Chambers et al. 2009) and the 
strong wing emission in the CO (3 -+ 2) line, indicating the presence of a powerful molecular 
outflow, all indicate star formation activity in the northernmost region of the filament. Due 
to its spatial location and mass, YSO number 29 in Shepherd et al. (2007) appears to be a 
good candidate for being the driver of the MM3-outflow. It is the most massive stellar object 
within MM3, with a mass of ~4.8-6.8 Mq and a luminosity of ~170-340 Lq. However, this 
is an intermediate mass star that may not be able to drive alone the massive outflow. In the 
same region there are three more stars of intermediate mass (less massive than the number 
29) with maximum masses between 4-5 Mq. Interferometric observations are needed to 
clarify if the MM3-outflow corresponds to an overlap of a few individual outflows, such as 
that seen in the central region of G34.43+0.24, or it is driven by a single star. 



5. Conclusions 

We carried out a multi-line study of cores within the fllamentary IRDC G34.43+0.24. 
The emission from the fllament was fully mapped in the CO (3 — ?■ 2) and ^^C0(3 — > 2) lines. 
In addition, emission from the central region of the fllament was observed in eight different 
molecular transitions [0^^0(3 ^ 2), CS(7 ^ 6), CS(2 ^ 1), C34S(2 -^ 1), HC0+(1 ^ 0), 
Hi3C0+(l -^ 0), SiO(2 -> 1) andCH30H(2 -^ 1)], and from the north region in the CS (2 -^ I] 
and 0^*^0(3 -+ 2) lines. The results of this study are summarized as follows. 
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The spatial distribution of the molecular emission in all, except methanol, lines is sim- 
ilar to that of the dust continuum emission derived from observations with similar angular 
resolution, showing a filamentary structure and four cores. 

The mass of the molecular cores derived assuming LTE conditions are 300 and 1460 Mq 
for the MMl and MM2 cores, respectively. The masses derived assuming virial equilibrium 
are 1100, 1500 and 1400 M© for the MMl, MM2 and MM3 cores, respectively. The average 
molecular hydrogen densities of the cores are 1.8 x 10^ cm~^ in MMl, 1.7x10^ cm"^ in MM2 
and 1.9 x 10^ cm~^ in MM3. We also find that the molecular hydrogen column densities at 
the peak positions of MMl and MM2 are, respectively, 2 x 10^^ cm"^ and 8 x 10^^ cm~^. 
The derived parameters of the MMl and MM2 cores are typical of massive and dense cores 
harboring high-mass YSO's (Plume et al. 1997; Garay et al. 2007). Those of the MM3 core 
suggest that it harbors young intermediate mass stars, in accord with the result obtained 
by Shepherd et al. (2007) from the modeling of the SED of the NIR sources detected using 
Spitzer. 

We modeled the molecular emission in the ^^00(3 -^ 2), C^^0(3 -^ 2), and CS(7 -^ 6) 
lines toward the MM2 core, concluding that this core is undergoing infalling motions, with 
an infall velocity of ~1.3 km s^^ and a mass infall rate of ~1.8 x 10^^ M© yr^^. This value 
is large enough to allow the formation of massive stars by accretion (see Osorio, Lizano & 
D'Alessio 1999). 

We report the discovery of a molecular outflow associated with the massive, dense MM3 
core located in the northern region of G34.43+0.24, and possibly of a second one associated 
with the MM4 core. The molecular outflow within MM3 has a total mass of 41 M©, a 
momentum of 270 M© km s~^, and a kinetic energy of 1120 Mq km^ s~^; indicating that 
it is driven by high-intermadiate mass stars. We also detected high velocity molecular gas 
toward the MMl and MM2 cores with total masses of 25 and 80 Mq, momenta of 230 and 
690 Mq km s~^, and kinetic energies of 1290 and 3630 Mq km^ s~^, respectively. 

The authors gratefully acknowledge support from CONICYT projects FONDAP No. 
15010003 and BASAL PFB-06. 
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Fig. 1. — Observed spectra of the ^^C0(3 — > 2) line emission toward G34.43+0.24, taken 
with the APEX telescope. The angular separation between panels is 20". Offsets are from 
the reference position at a = 18''53™17."40 and 6 = 01°24'55" (J2000). In each box the 
velocity scale ranges from 40 to 80 km s~^. The antenna temperature scale is from -2 to 12 
K. 
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Fig. 2. — Contour maps of velocity-integrated ambient molecular emission toward the central region of 
G34.43+0.24, overlaid on a 1.2 mm dust continuum emission map (gray scale). The range of integration in 
visr is from 52.9 to 61.1 km s^^. The FWHM beam is shown in the bottom left corner. Contour levels are 
from 20% to 90% (in steps of 10%) of the peak integrated intensity of 49.9 K km s^^ for "C0(3 -^ 2) (a); 
7.7 K km s-i for CS(7 -^ 6) (&); 16.9 K km s^^ for CS(2 -^ 1) (c); 3.4 K km s^^ for SiO(2 -^ 1) (d); 3.7 K 
km s-i for Hi3C0+(l -^ 0) (e); and 14.3 K km s^^ for Ci**0(3 -> 2) (/). The blue lines indicate the 50% 
contour level. The gray scale levels are 1, 2, 3, 5, 7, 10, 15, 20 and 30 x 0.12 Jy beam~^. The triangles 
mark the peak positions of the MMl, MM2 and MM4 dust cores. 
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Fig. 3. — Contour maps of velocity-integrated ambient CH3OH emission toward the central region of 
G34.43-I-0.24, overlaid on a 1.2 mm dust continuum emission map (gray scale). The FWHM beam is shown 
in the bottom left corner. Contour levels are from 20% to 90% (in steps of 10%) of the peak integrated 
intensity of 5.1 K km a~^ for CH30H(2o -^ lo) A+ (a); 4.1 K km s^^ for CH30H(2_i -^ l_i) E (6); and 
1.1 K km s^^ for CH30H(2o — s> lo) E (c). The blue lines indicate the 50% contour level. The gray scale 
levels are 1, 2, 3, 5, 7, 10, 15, 20 and 30 x 0.12 Jy beam^^. The triangles mark the peak positions of the 
MMl, MM2 and MM4 dust cores. 



-26- 



r29'00" - 



r28'00" - 



r27'00" r 



o r29'00' 

o 

o 

c\l 



c 

o 

D 

0) 
Q 



18'00' 



'00' 



r29'00" - 



1°28'00" - 



r27'00" r 




h m s 

18 53 20 



15^ 



Right Ascension (J2000;; 



Fig. 4. — Contour maps of velocity-integrated molecular line emission toward the MM3 core (Garay et 
al. 2004), overlaid on a 1.2 mm dust continuum emission map (gray scale). The range of integration in visr 
is from 56.1 to 61.5 km s^^. Contour levels are from 20% to 90% (in steps of 10%) of the peak integrated 
intensity of 12.5 K km s^^ for 13C0(3 ^ 2) (a) and 3.7 K km s^^ for CS(2 ^ 1) (6); and from 3 to 21 ct 
in steps 3 a, with a = 0.14 K km s-\ for Ci*0(2 -^ 1) (c). The gray scale levels are 1, 2, 3, 5, 7, 10, 15, 
20 and 30 x 0.12 Jy beam^^. The blue lines indicate the 50% contour level. The triangle marks the peak 
position of the MM3 core. The FWHM beam is shown in the bottom left corner. 
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Fig. 5. — Molecular spectra observed toward the peak position of the millimeter continuum 
of MMl core. Transitions and telescopes are given in the top corner of the spectra. The 
vertical dotted line indicates the systemic velocity of the ambient gas of 57.6 km s~^ (Miralles, 
Rodriguez & Scalise 1994). 
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Fig. 6. — Molecular spectra observed toward the peak position of the millimeter continuum 
of MM2 core (IRAS 18507+0121). Transitions and telescopes are given in the top corner of 
the spectra. The vertical dotted line indicates the systemic velocity of the ambient gas of 
57.6 km s~^ (Miralles, Rodriguez & Scalise 1994). 
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Fig. 7. — Molecular spectra observed toward the peak position of the millimeter continuum 
of MM3 core. Transitions and telescopes are given in the top right corner of the spectra. 
The vertical dotted line indicates the systemic velocity of the ambient gas of 58.7 km s"^ 
(estimated from the optically thin i3C0(3 -^ 2), 0^^0(2 ^ 1) and CS(2 -^ 1) lines). 
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Fig. 8. — Position-velocity diagrams along lines in the north-south direction (i.e., parallel to the filament). 
Transitions and telescopes are indicated in the top left corner of each panel. Left: Cut passing through MMl 
or position (0", 20"). Contour levels are 3, 9, 15, 23, 31, 41, 51, and 61 x 0.27 K for C0(3 -^ 2); 3, 6, 9, 
12, 17, 22, 28, and 34 x 0.08 K for CS(2 -^ 1); 3, 4, 6, and 8 x 0.07 K for SiO(2 -^ 1); and 1.5, 3, 6, 9, 
12, 15, and 19 x 0.37 K for 13C0(3 -^ 2). Right: Cut passing through MM2 or position (20", 0"). Contour 
levels are 3, 8, 13, 18, 28, 38, 48, and 58 x 0.27 K for C0(3 -^ 2); 3, 6, 11, 16, 21, 28, 35, and 42 x 0.08 K 
for CS(2 -^ 1); 3, 5, 7, and 9 x 0.07 K for SiO(2 -^ 1); and 1.5, 3, 6, 10, 14, 19, 25, and 31 x 0.37 K for 
i3CO(3 ^ 2). 
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Fig. 9. — Contour maps of velocity- integrated wing emission toward the central region of G34.43-I-0.24, 
in different molecular lines. Transitions and telescopes are indicated in the top left corner of each panel. 
The arrows in the panels indicate where are the outflows corresponding to each core. Contour levels for the 
blueshifted emission are 20, 35, 50, 80, and 110 a for C0(3 -^ 2); 3, 5, 7, 9, 12, and 16 a for SiO(2 -^ 1); 
and 3, 9, 15, 25, and 40 a for CS(2 -^ 1). For the redshifted emission, these are 20, 35, 50, 80, 110, 150, 190, 
and 230 a for C0(3 -^ 2); 3, 5, 7, 9, 12, 16, 24, and 32 a for SiO(2 -^ 1); and 3, 9, 15, 25, 40, 60, 85, and 
110 a for CS(2 -^ 1). The a values are: 0.30 K km s^^ for the CO line, and 0.04 K km s^^ for the SiO and 
CS lines. 
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Fig. 10. — Position- velocity diagrams of CO (3 -^ 2) and ^^C0(3 — )► 2) toward the north- 
ernmost region of G34. 43-1-0.24 along a line running since the north-south direction passing 
through of MM3 (position (40", 0") in Figure [1]). Transitions and telescopes are indicated 
in the top left corner of each panel. The contour levels for C0(3 — )■ 2) are 3, 6, 10, 15, 22, 
29 and 36 a, with a = 0.27 K. For 1^00(3 -^ 2) are 1.5, 3, 6, 8 and 10 a, with a = 0.37 K. 
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Fig. 11. — Contour maps of velocity- integrated wing emission of CO (3 — ?■ 2) and 
^^C0(3 — > 2) toward the northernmost region of G34.43+0.24. Sohd hnes represent 
blueshifted emission and dashed hnes represent redshifted emission. The asterisk shows 
the position of the MM3 core. Top: CO (3 — > 2) emission. Contour levels for the blueshifted 
emission are 20, 30, 40, 50, 75 and 100 a {a = 0.35 K km s~^); and for redshifted emission, 
these are 20, 30, 60, 90, 120 and 150 a {a = 0.25 K km s"^). Bottom: ^^00(3 -^ 2) emission. 
Contour levels for the blueshifted emission are 3, 4, 5, 7 and 9 a (a = 0.23 K km s~^); and 



for redshifted emission, these are 3, 4, 5, 6 and 7 a {a 



0.19 K km s-i). 
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Fig. 12.— ^3Qo(3 -> 2), 0^^0(3 ^ 2), and CS(7 -^ 6) spectra of the MM2 pealc position 
with the simple two layer model fits of Myers et al. (1996). The solid lines show the observed 
spectra and the dashed (red) lines represent the model. The values derived by the model are 
shown in the right side of the figure. 
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Table 1: Properties of 1.2 mm Cores 





Peak Position 


FWHM Diameter 
Angular Physical 






a(J2000) 


5(J2000) 


Mass 


Core 






(") 


(pc) 


M© 


MMl 


18 53 18.0 


01 25 24 


16 


0.19 


1187 


MM2 


18 53 18.6 


01 24 40 


26 


0.42 


1284 


MM3 


18 53 20.4 


01 28 23 


24 


0.38 


301 


MM4 


18 53 19.0 


01 24 08 


24 


0.38 


253 



Note. — Units of right ascension are hours, minutes and seconds, and units of declination are degrees, 
arcniinutes and arcseconds. 

Note. — Properties adopted from Rathborne, Jackson & Simon (2006). 
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Table 2: Summary of Observational Parameters 





Frequency 


Beam (FWHM) 


Spacing 


Observed 


Region Size 


Av 


rms Noise (Ta) 


Observation 


(GHz) 


(aresee) 


(arcsec) 


Positions 


(arcsec^) 


(km s-i) 


(K) 


APEX 


i»CO(3^2)... 


330.588 


18.4 


20 




145=^ 


120 X 700 


0.11 


0.37 


CO(3 ^2) 


345.796 


17.6 


20 




145^^ 


120 X 700 


0.11 


0.27 


CS{7 ^6) 


342.883 


17.8 


20 




21^ 


60 X 140 


0.11 


0.18 


Ci»0(3^2)... 


329.331 


18.5 


20 




30'' 


100 X 120 


0.11 


0.29 


Nobcyama 


CS(2 ^1) 


97.981 


15.5 


20 




30'' 


100 X 160 


0.11 


0.08 


SiO(2^1) 


86.846 


17.7 


20 




30^ 


100 X 160 


0.13 


0.07 


C34S(2^1) ... 


96.412 


15.8 


20 




SI*" 


100 X 160 


0.12 


0.10 


HCO+(l ^0).. 


89.188 


17.3 


20 




31'' 


100 X 160 


0.12 


0.08 


Hi3cO+(l->0) 


86.754 


17.7 


20 




30^ 


100 X 160 


0.13 


0.07 


CH30H(2 ^ 1) 


96.744 


15.7 


20 




31*' 


100 X 160 


0.12 


0.09 


SEST 


CS(2->1) 


97.981 


51.0 


30 




25<= 


100 X 100 


0.13 


0.07 


Ci80(2^- 1)... 


219.560 


24.1 


30 




25<= 


100 X 100 


0.12 


0.18 



"Covering the whole filamentary structure. 

''Covering the central region of G34.43+0.24, including the MMl, MM2 and MM4 cores. 

^Covering the MM3 core, in the northern region of the G34.43+0.24 filament. 
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Table 3. Observed Parameters of Molecular Cores 



Molecular Angular Size 

Core Transition (" x ") 

MMl CS(7^6) 37.2x20.6 



Vlsr 


AV 


/ Tmbdv 


Physical Size 


(km s-i) 


(kms-i) 


(K km s-i) 


(pc2) 


57.38 ± 0.21 


5.69 ± 0.71 


6.94 ± 0.80 


0.66 X 0.38 



CS(2^1) 45.4x33.0 58.13 ± 0.01 5.29 ± 0.02 15.10 ± 0.04 0.82x0.60 

Hi3C0+(l -^ 0) 33.0 X 24.8 58.36 ± 0.02 2.75 ± 0.06 2.78 ± 0.05 0.60 x 0.44 

MM2 CS(7^6) 37.2 x 24.8 57.96 ± 0.03 5.27 ± 0.09 9.41 ± 0.13 0.66 x 0.44 

CS(2^1) 41.2x37.2 57.35 ± 0.01 5.19 ± 0.03 15.41 ± 0.07 0.74x0.66 

H"CO+(1^0) 41.2x33.0 57.88 ± 0.02 3.96 ± 0.05 3.52 ± 0.04 0.74x0.60 

MM3 "Qo(3^2)... 74.2 x 53.6 58.77 ± 0.02 3.01 ± 0.04 10.31 ± 0.11 1.34 x 0.96 

CS(2^1) 86.6x70.0 58.73 ± 0.01 3.43 ± 0.03 3.44 ± 0.02 1.54x1.24 
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Table 4: Derived LTE Core Parameters 



Core 


Transition 


T-j. 


(cm-2) 


N(H2) 
(cm-2) 


Mlte 

(Mo) 


MMl 


Ci80(3 ^ 2) 
i3CO(3 ^ 2) 


0.10 
0.75 


5.2 X 10^5 
3.9 X 10^6 


2.0 X 10^2 


330 


MM2 


Ci80(3 ^ 2) 
i3CO(3 ^ 2) 


0.18 
1.39 


1.1 X 10^6 
8.7 X 10^6 


4.3 X 10^2 


1460 















Note. — Velocity range for the MMl core: 54.6-59.3 km s ^. 
Velocity range for the MM2 core: 53.4-59.3 km s^^. 
° Column density for each molecule. 
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Table 5: Derived Virial Core Parameters 



Radius^ AV M^^^ n(H2) 

Core Transition (pc) (km s~^) (Mq) (cm~^) 

MMl CS(7^6) (125 5^69 1700 3.9 x 10^ 

CS(2-^1) 0.35 5.29 2060 1.7 x 10^ 

Hi3C0+(l -^ 0) 0.26 2.75 410 8.3 x 10^ 

MM2 CS(7^6) (127 5^27 1580 2.8 x 10^ 

CS(2^1) 0.35 5.19 1980 1.6 x 10^ 

Hi3co+(1^0) 0.33 3.96 1090 1.1x10^ 

MM3 i3CO(3^2)... 057 3M 1090 2.1 x 10^ 

CS(2^1) 0.69 3.43 1710 1.8x10^ 



''Geometric mean of the observed semi-major and semi-minor axes. 
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Table 6: Core Masses from Different Authors 





LTE Mass 
(This work) 

(2) 


Virial Mass 

(This work) 

(3) 


LTE Mass 
Shepherd et ah (2007) 

(4) 


Dust Mass (1.2 mm) 


Core 

(1) 


Garay et at (2004) Rathborne et al. (2006) 
(5) (6) 


MMl 
MM2 
MM3 


330 
1460 


1130 
1510 
1370 


75 
690 


550 1187 
2700 1284 
780 301 












Note. 


— Masses are 


in Mp,. 
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Table 7: Core Parameters at Peak Positions 









N^ 


N(H2) 


Core 


Transition 


Tu 


(cm-2) 


(cm-2) 


MMl 


Ci80(3 ^ 2) 


0.04 


2.3 X 10^5 


8.5 X 10^1 




i3CO(3 -^ 2) 


0.30 


1.7 X 10^6 






C34S(2 ^ 1) 


0.02 


1.9 X 10^3 


4.3 X 10^2 




CS(2 ^ 1) 


0.41 


4.3 X 10^^ 




MM2 


CiS0(3 ^ 2) 


0.37 


2.3 X 10^6 


8.5 X 10^2 




i3CO(3 ^ 2) 


2.80 


1.7 X 10^^ 






C34S(2 ^ 1) 


0.03 


3.5 X 10^3 


8.0 X 10^2 




CS(2 ^ 1) 


0.69 


8.0 X 10^^ 





Note. — Velocity range for the MMl core: 54.0-59.3 km s ^. 
Velocity range for the MM2 core: 53.5-59.2 km s^^. 
° Column density for each molecule. 
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Table 8. Derived Parameters of Molecular Outflows 



Core Lobe 



Mlobe Mhiddcn 

(Mo) (Mo) 



Mtotal 

(Mo) 



(Mo km s 



-l^ 



"pmin 



"pmax 



(Mo km^ s 



2 c-2\ 



MMl 


Blueshifted 


3 


3 


6 


2.8 X 10^ 


1.1 X 10^ 


9.2 X 10^ 


9.2 X 10^ 




Redshifted 


13 


6 


19 


9.5 X 10^ 


3.3 X 10^ 


3.5 X 102 


2.9 X 10^ 




Total 


16 


9 


25 


1.2 X 10^ 


4.4 X 10^ 


4.4 X 10^ 


3.8 X 10^ 


MM2 


Blueshifted 


15 


21 


36 


1.8 X 10^ 


5.5 X 10^ 


6.4 X 102 


4.2 X 103 




Redshifted 


32 


12 


44 


2.2 X 10^ 


6.7 X 102 


7.8 X 10^ 


5.2 X 103 




Total 


47 


33 


80 


4.0 X 10^ 


1.2 X 10^ 


1.4 X 10^ 


9.4 X 10^ 



MM3 Blueshifted 10 
Redshifted 8 

Total 18 



14 


24 


6.8 X 10^ 


3.5 X 10^ 


1.8 X 10^ 


2.5 X 10^ 


9 


17 


4.9 X 10^ 


2.4 X 10^ 


1.1 X 102 


1.8 X 10^ 


23 


41 


1.2 X 10^ 


5.9 X 10^ 


2.9 X 10^ 


4.3 X 10^ 



